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This  report  describes  the  results  of  a  two-dimensional,  thick, 
turbulent  boundary^ layer  separation  study  conducted  at  the  Naval 
Ordnance  Laboratory.  The  purpose  of  this  investigation  was  to 
obtain  pressure  distribution  data  in  the  separation  region  ahead 
of  the  forward  facing  step  '{used  to  separate  the  boundary  layer) 
for  the  case  where  the  step  height  was  less  tKan.  the  .boundary^-layer 
•thickness.  This  study  was  performed'  for  the  Naval  Air  Systems 
Command  under  Task  Number  A320  320C/004B/WP32-322-202 . 
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INTRODUCTION 

This  report  describes  the  experimental  results  of  a  thick,  two- 
dimensional  turbulent  boundary-layer  separation  study  conducted  at 
the  Naval  Ordnance  Laboratory  (NOL)  .  The  separation  was 
induced  by  a  forward  facing  step.  The  purpose  of  this  investigation 
was  to  obtain  pressure  distribution  data  in  the  separation  region 
ahead  of  the  step  for  the  case  where  the  step  height  was  less  than 
the  boundary- layer  thickness.  This  information  is  technically 
important  to  those  studying  boundary- layer  separation  and-  has  a 
variety  of  possible  applications  which  include  supersonic  flight 
controls.  Numerous  investigators  (Refs.  (1)',  (2)  ,  (3)  and  (4) ) 
have  studied  separation  phenomena  ahead  of  forward  facing  steps 
for  the  turbulent  boundary-layer  case.  Of  these.,  the  work  of 
Zukoski  (Ref.  (1))  is  particularily  good  since  it  is  both  a 
review  and  an  analysis  of  numerous  two-dimensional  experimental 
works.  As  a  result  of  these  reported  works,  simple  correlations 
have  been  determined  for  the  upstream  pressure  distributions  and 
for  the  s-'  ^e-force  calculations  of  the  forward  facing  step  model 
limited  to  the  case  for  h/6  >  1.5. 

In  a  recent  jet  interaction  control  s^ady  at  NOL,  (Ref.  (5)), 
the  similarity  of  the  upstream  pressure  distributions  for  both 
the  step  and  interacting  jet  separation  regions  was  demonstrated. 
This  similarity  established  the  fact  that  jet  and  step  separation 
phenomena  are  analogous.  The  study  of  Reference  (5)  was  limited 
to  the  case  where  the  boundary  layers  were  thin  relative  to  the 
height  of  the  particular  obstruction  which  induces the  boundary - 
layer  separation.  Hence,  the  thick  boundary- layer  case  (h  <  6) 
remained  to  be  investigated  and  is  the  subject  of  this  report. 

SYMBOLS 

D  see  Equation  (4),  inch 

F  side  force,  Ibf/inch 

g  any  dimensional  flow-field  property 

G  any  dimensionless  flow-field  parameter 

h  step  height,  inch 

M  Mach  number 

n  exponent  of  the  turbulent  boundary-layer  profile 

P  pressure,  psia 

PI  first  plateau  pressure,  psia 

R  gas  constant,  ft-lbf/lbm  °R 

Re^  Reynolds  number  based  upon  boundary- layer  thickness 

T  temperature,  °R 

X  distance,  measured  along  the  flat  plate  surface,  inch 
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SYMBOLS  .{Continued) 

projected  shock,  distance/  see  sketch  Page  3./  inch 

separation,  distance,  see  sketch  Page  8,  inch 

initial  steep  pressure  rise  distance,  see  sketch  Page  8 
inch 

see  Figure  14d' 
see .Figure  14d 
specific  heat  ratio, 
boundary- layer  thickness,,  inch 
dynamic  viscosity ,  slug/ft  sec 
density,;  lom/ft3 


CO 

remote  free  stream 

0 

stagnation 

oil 

oil  flow  study 

EXPERIMENTAL  APPARATUS 

The  experimental  tests  were  performed  in  the  NOL  Boundary 
Layer  Channel  (Ref.  (6)).  A  schematic  of  the  Boundary  Layer  Channel 
is  presented  in  Figure  1.  The  flexible  nozzle  plate  is  directly 
opposed  to  a  flat  test  plate  along  which  thick  boundary  layers  are 
generated  for  testing.  The  flat  test  plate  is  106  inches  long 
and  14  inches  wide.  The  wind-tunnel  nozzle  side  walls  are 
hinged  to  the  flat  test  plate  and  serve  as  access  doors  to  the 
test  section.  In  the  closed  position,  the  doors  are  sealed  by 
pressing  against  sealing  gaskets  located  in  the  edges  of  the  flat 
and  flexible  plates.  Glass  viewing  ports  were  located  in  the  doors 
thereby  permitting  one  to  view  the  flow  field. 

The  flexible  nozzle  plate  is  held  in  position  by  13  screw 
jacks.  The  particular  contour  selected  for  the  flexible  plate 
determines  the  pressure,  gradient  along  the  length  of  the  wind- 
tunnel  nozzle.  These  gradients,  of  course,  directly  influence  che 
test  boundary  layer  generated  along  the  flat  plate.  For  this 
experimental  study,  a  constant  pressure  gradient  contour  equal  to 
zero  was  specified  for  the  test  section.  This  contour  was 
computed  using  the  method  of  characteristics  which  included 
corrections  for  the  boundary- layer  displacement  thickness. 

Figure  2a  presents  the  experimental  Mach  numbers  obtained  in  the 
test  rhombus  for  the  various  wind-tunnel  supply  pressures  while 
Figure  2b  presents  the  boundary-layer  thickness  distribution  along 
the  plate.  The  Mach  number  fluctuations  in  Figure  2a  ranged  from 
2.5  percent  to  0.5  percent  as  the  supply  pressure  varied  fro*”  1  to  5 
atmospheres.  The  overall  change  of  the  Mach  number  in  th*s 
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pressure  range  was  about  4.8  percent.  This  overall  Mach  number 
dhange  was  attributed  to  changes  in  the  boundary-layer  displacement 
thickness.  The  local  fluctuations  in  Mach  numbers  were  attributed 
to  measurement  errors  with  the  maximum  errors  occurring  at  the 
lowest  stagnation  pressures. 

Figure  3  is  an  oblique  view  of  the  step  model  and  the 
installed  glass-ported  end  plates.  These  end  plates  (details 
are  given  in  Fig.  4)  were  mounted  eight  inches  apart,  symmetrically 
about  the  centerplane  of  the  half  nozzle  at  zero  angle  of  attack. 

The  forward  most  tips  of  the  two  end  plates  were  located  at 
the  76.0-inch  station  (measured  from  the  sonic  throat  station  of 
the  wind-tunnel  nozzle) „  This  distance  aligned  the  various 
glass  windows  properly  for  viewing.  The  metal  retaining 
ring  which  held  the  door's  glass  window  in  place,  blocked  the 
flow-field  view  for  the  first  0.25  inch  off  the  test  plate 
surface  as  seen  in  Figure  4 . 

A  simple,  Z  configuration,  spark  shadowgraph  system  was 
used  to  record  all  photographic  data. 

The  primary  experimental  data  were  the  upstream  pressure 
distributions.  Figure  5  shows  the  physical  arrangement  of  the 
static  pressure  taps.  The  numbers  shown  on  the  right-hand  side 
of  the  taps  are  the  hole  identification  numbers;  those  on  the 
left-hand  side  -re  the  hole  distances  measured  from  the  step 
mounting  plar-.  .  The  two  test  locations  of  the  step  face  were 
located  0.375  1  0.6875  inches  ahead  of  this  step  mounting  plane. 

A  total  of  wcren  steps  were  tested.  The  installation  detail 
and  a  listing  of  the  step  heights  tested  are  given  in  Figure  6. 

To  obviaue  any  possibility  that  the  high-pressure  gas  found  in 
the  forward  separation  bubble  could  leak  out,  black  plastic 
tape  w7as  used  to  seal  all  interior  cracks  and  corners.  These 
seals  are  visible  in  the  photograph  of  Figure  3. 

To  measure  the  static  pressure  distributions,  eight  strain- 
gage  pressure  transducers  were  directly  connected  to  eight 
scanner  valves.  These  scanner  valves  were  switched  sequentially 
through  12  valve  positions  according  to  the  pressure  tap  array 
of  Table  1.  Note  that  the  valve  position  1  was  reserved  for 
calibration  purposes.  This  array  identifies  the  various  pressure 
tap  locations  (Fig.  5)  with  the  corresponding  scanner  valve  and 
sequence  position.  Tap  numbers  65  through  78  were  not  used 
for  this  step  study.  The  electrical  switching  equipment  was 
designed  to  permit  only  forward  sequential  switching  of  the 
scanner  valves.  The  switch  position  was  sensed  and  recorded  along 
with  each  sample  of  pressure  data  thus  eliminating  error  with  regard 
tc  the  static  pressure  tap  matrix.  Steady  tunnel  conditions 
were  maintained  over  the  two-to-four-minute  time  period  required 
to  measure  all  static  pressures.  Instantaneous  wind-tunnel 
conditions  were  used  for  all  data  reduction. 
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TABLE  1 

STATIC  PRESSURE  TAP  IDENTIFICATION  ARRAY 


SCANNER  VALVE  SEQUENCE  POSITION 


The  transducer  signals  were  recorded  using  NOL's  updated 
"Portable  Automatic  Data  Recording  Equipment,"  (PADRE),  an  earlier 
version  of  which  is  described  in  Reference  (7) .  This  recording 
unit,  shown  in  Figure  7,  provides  a  visual  display  of  the  data 
output.  These  data  are  recorded  on  IBM  cards  which  are  then  used 
in  an  IBM  7090  computer  for  the  final  data  reduction. 

The  final  piece  of  equipment  used  in  this  study  was  the 
boundary-layer  probe  (See  Appendix  B  for  further  information  and 
reference  works) .  This  boundary- layer  survey  mechanism  was 
externally  attached  to  the  flat  test  plate  of  the  wind-tunnel 
nozzle  (shown  in  Fig.  7)  at  each  of  two  available  probing  ports. 
These  probing  ports  were  located  on  the  wind-tunnel  nozzle's 
centerplane,  48.0  and  60.0  inches  downstream  of  the  wind-tunnel 
throat.  The  boundary-layer  description  as  determined  from  the 
probe  measurements  is  presented  in  Appendix  B. 

DIMENSIONAL  ANALYSIS 

The  test  program  was  organized  using  the  technique  of 
dimensional  analysis.  Figure  8  is  a  schematic  representation  of 
the  wind-tunnel  nozzle-step  geometry.  'or  this  analysis  it  is 
assumed  that  no  heat-transfer  or  chemical  reactions  take  place. 
Consideration  of  the  associated  wind-tunnel  variables  and  the 
simple  two-dimensional  flat  plate-step  geometry  results  in  the 
following  general  dimensional  relationship  for  any  flow-field 
property,  g. 


g  =  g  [u,  P,  p,  R,  y,  ji,  <5,  h) 


(1) 
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An  analysis  of  the  variables  and  primary  dimensions  in  Equation  (1) 
above  reveals  that  a  nondimens ional  parameter,  G,  may  be  determined 
which  will  be  a  function  of  at  least  four  dimensionless  numbers.. 

This  can  be  written  as  ' 

G  =  G  [Reg,  M,  y,  h/61  (2) 

Since  both  y  =  1.4  and  M  =  4.9  were  constants  for  this  test. 

Equation  (2)  above  reduces  to 

G  =  G  [Reg ,  h/S]  (3) 

Equation  (3)  expresses  the  fact  that  ly  dimensionless  flow-field 
property,  G,  (as  for  instance  the  upscream  pressure  signature 
">f  the  separated  flow  region)  is  to  be  investigated  as  a  function 
of  the  flow  geometry,  h/6,  and  the  Reynolds  number  based  upon  the 
boundary- layer  thickness,  Reg.  As  seen  in  Equation  (3),  the 

significant  length  to  be  associated  with  the  fluid  was  the 
bo-ndary-layer  thickness,  6,  measured  at  the  point  where  the 
boundary- layer  separates  from  the  flat  plate.  Normally,  one 
may  also  choose  as  fluid  lengths  either  the  flat  plate  distance 
measured  from  a  leading  edge  or  the  distance  measured  from  the 
end  of  transition  (turbulent  case).  However,  for  this  test, 
the  first  length  did  not  exist  (by  wind-tunnel  design) ,  and 
the  second  length  was  unknown.  Even  if  the  precise  location  of 
the  end  of  transition  had  been  known,  though,  the  significance 
of  the  distance  measured  from  this  point  would  be  unclear  as 
transition  certainly  occurred  within  the  accelerating  flow  region 
of  the  wind-tunnel  nozzle  and  not  in  the  test  rhombus. 

EXPERIMENTAL  PROGRAM  AND  PROCEDURE 

The  experimental  test  program  outline  is  given  in  Table  2. 

This  outline  follows  the  functional  expression  given  by  Equation  (3) . 

TABLE  2 


PRESSURE  TEST  PROGRAM* 


foot.  Included  in  this  study  was  a  0.130-inch  high  step  not  listed 
above . 
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The  Boundary  Layer  Channel  contour  was  checked  prior  to 
testing  fay  experimentally  measuring  the  pressure  gradient  throughout 
the  test  rhombus..  Several  preliminary  pressure  surveys  and  contour 
adjustments  were  necessary  to  produce  satisfactory  results. 
Subsequently,  the  glass-ported  end  plates  were  installed  arid  the 
final  pressure  gradients  were  measured  using  the  surface  static 
pressure  taps  located  between  the  glass  end  plates.  The  results 
of  this  survey  are  presented  in  Figure  9. 

The  procedure  followed  to  record  the  pressure  gradient  da.ta 
was  identical  with  that  followed  later  to  record  the  pressure 
distributions  ahead  of  the  step.  Quite  simply,  this  procedure 
required  a  sequential  shifting  of  the  scanner  valves,  a  pauoc 
to  observe  steadiness  of  the  pressure  transducer  outputs,  and 
a  recording  operation.  One  can  identify  the  time  sequence 
followed  for  taking  data  and  also  identify  the  particular  pressure 
tap-transducer  combination  used  by  comparison  of  Table  1  with 
Figure  5. 

Two  pressure  distributions  were  recorded  for  each  step,  one 
distribution  for  each  of  the  two  available  step  mounting  positions 
(see  Fig.  6) .  This  provided  offset  pressure  distributions  for 
comparison  purposes  and  assured  repeatibility  since  separate 
wind-tunnel  startups  were  involved.  Figure  10  shows  a  typical 
example  of  offset  pressure  distributions.  Excellent  repeatibility 
was  obtained  for  all  cases.  Hence,  only  one  distribution  is 
presented  in  this  report. 

EXPERIMENTAL  RESULTS 

The  significant  experimental  test  parameters  determined  by 
dimensional  analysis  were  Re^  and  h/6.  However,  the  more 

practical  dimensional  parameters,  Re/L  and  h,  were  used  experimentally 
as  the  independent  test  variables.  No  problem  was  cheated  with 
this  simplification  since  <5  varied  only  fractionally  over  the 
range  of  Reynolds  numbers  covered  in  this  experiment,  (Fig.  2b) . 

Hence,  variations  in  Re^  follow  almost  directly  the  experimenta.l 

variations  in  Re/L  and  similarly  so  for  h/S  and  h. 

UPSTREAM  PRESSURE  DISTRIBUTION 

Figure  11  presents  three  sets  of  pressure  distributions. 

The  data  shown  in  Figures  11a  and  11b  were  taken  at  Reynolds 
numbers  of  4.0  x  10 e  per  foot  and  0.8  x  10s  per  foot,  respectively, 
over  the  test  range  of  step  heights.  In  contrast,  Fiaure  11c 
presents  data  taken  at  the  constant  step  height  of  1.400  inches 
over  the  experimental  range  of  Reynolds  numbers  (Table  1) .  As 
indicated  by  the  data,  the  step  height  influences  both  the  extent 
of  separation  and  the  value  of  the  maximum  pressure  at  the  first 
pressure  plateau,  whereas  the  Reynolds  number  per  foot  is  observed 
to  influence  the  pressure  profile  in  the  upstream  region  near  the 
beginning  of  the  pressure  rise. 
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The  small  amount  of  scatter  observed  in  the  data  in  Figure  lib 
may  be  attributed  to.  the  very  low  values  of  static  pressures  which 
existed  in  the  flow  field  for  those  test  conditions.  Pressures  as 
low  as  0.03  psia  were  measured  using  0-5.0  psia  full-scale  pressure 
transducers  calibrated  over  a  reduced  pressure  range  of  0-1.0  psia. 
Thus,  the  pressures  of  0.03  ps..a  were  measured  near  the  practical 
lower  limit  of  the  pressure  transducer  and  some  fluctuations  were 
inevitable . 

OIL  FLOW  RESULTS 

Figure  12  presents  all  of  the  oil  flow  photographic  data 
obtained.  Silicon  oil  {300  centistokes) ,  with  white  Titanium 
oxide  powder  in  suspension,  was  painted  in  vertical  strips  on  the 
flat  test  plate  of  the  Boundary  Layer  Channel  (the  wind-tunnel 
nozzle  is  oriented  vertically) .  After  operating  the  wind  tunnel 
for  a  few  minutes  the  oil  flow  pattern  was  developed.  When  the 
wind  tunnel  was  shut  down,  the  accumulation  of  oil  located  at  the 
separation  line  flowed  (under  the  influence  of  gravity)  downward 
on  the  flat  test  plate.  This  explains  the  rather- wide  separation 
lines  seen  in  Figure  12.  Separation  distance  was  measured  from  the 
step  face  to  the  topmost  edge  of  this  line.  One  rather  interesting 
finding  of  this  oil  flow  study  was  the  result  of  Figure  13.  By 
making  all  measured  oil  flow  points  coincident,  it  appears  that 
a  single  pressure  signature  exists  with  local,  individual 
deviations  occurring  only  near  the  region  where  the  reattachment 
pressure  rise  exists.  Note,  however,  that  the  plateau  pressures 
are  also  different  for  each  step.  In  all  cases,  oil  flow 
separation  was  located  at  the  point  in  the  pressure  distribution 
where  P/P^  =  1.2.  This  compares  favorably  with  the  recently 
reported  results  of  P/P^  =  1.2  to  1.25  (Ref.  (8)). 

SELECTION  AND  CORRELATIONS  OF  THE  SCALING  LENGTHS 

In  Reference  (5)  pressure  distributions  ahead  of  steps  and 
jets  for  h/<$  >  1.5  were  obtained  by  plotting  (P-P^) / (Pl-P^)  versus 

X/XS.  In  these  correlations  the  plateau  pressure  rise,  (Pl-P^) ,  and 

the  known  separation  distance,  XS,  were  used  as  the  normalizing 
parameters.  Because  of  this  previous  success  identical  parameters 
were  selected  for  correlating  the  pressure  data  of  the  thick 
boundary-layer  case.  The  plateau  pressure  rise  data  were  taken 
from  the  pressure  distributions  of  Figure  11  and  the  separation 
distances  were  available  from  the  previously  mentioned  oil  flow 
study.  However,  the  use  of  these  normalizing  parameters  did  not 
produce  a  universal  pressure  profile.  Father,  the  diverging  nature 
of  the  plotted  results  indicated  that  rhe  fault  was  related  to 
the  measured  oil  flow  separation  distance.  Since  this  distance 
did  not  scale  the  pressure  distributions,  a  separation  criterion 
established  by  Zukoski  (Ref.  (1))  was  selected  for  use  hopefully 
to  obtain  a  scaling  length  which  would  correlate  the  data. 

According  to  this  criterion,  the  separation  pressure  rise  is  73 
percent  of  the  plateau  pressure  rise.  Application  of  this  ru±e  to 


N0LTR  7,3-98  • 

tharia?heSoiinis  uSfcnS:  •l0eatl°n  °f  *ST  m  means  other 

on>-^S^bSeq^Jnt  plots  using  these  data  established  a  universal 
vSs  of X  > -xf  n?h  °  *  X/XS  4  when  exteiSS  t6 

£  x  XS'  the  pressure  profile  again  diverged  into  a 
impliedff Tv!65’  H|ncf'  f  change  in  relevant  scale  length  is 
JSiJe  11  dll  The  P^ssure  distributions  of 

significant  to  £rof  which  an  appropriate  distance, 

t°  thG  *nitxal  steep  pressure  rise  region,  could  be 
earned  and  measured.  This  distance  was  designated  as  the  "ini-Hai 
steep  pressure  rise  distance" ^XS.  ^rgnarea  as  the  initial 

to  obtainfAXS°Win9  SkStCh  presents  the  graphical  method  used 


DISTANCE 


M 


NOLTR  73-98 


4'J 


ft 


* 


A 


V 


The  use  of  AXS  as  a  scaling  length  completed  the  universal 
pressure  distribution.  For  convenience,  in  the  presentation  of 
the  universal  pressure  distribution  curve  of  Figure  14,,  the 
separation  point  determined  by  Zukoski’s  criterion  (Ref.  (1))  as 
established  as  the  point  X  =  0,  permitting  both  sets  of  normaliz.ed 
distances  to  vary  from  0  to  1.0.  The  experimental  results  of  this 
correlation  are  presented  in  Figures  14a,  14b,  and  14c.  These 
three  distributions  appear  to  be  of  the  same  parent . population 
and  are  adequately  represented  by  the  curve  in  Figure  14d,  the 
Composite  Universal  Pressure  Distribution. 

The  existence  of  a  universal  pressure  profile  implies  that  the 
functional  dependence  of  the  pressure  distributions  with  Re/L  and 
h/6  (see  Fig.  11)  is  incorporated  in  the  scale  parameters  AXS  and 
XS. 


To  investigate  this  functional  dependence,  the  projected  shock 
distance,  XP  (see  previous  sketch),  was  introduced. 

As  seen,  XP  is  determined  by  the  straight  line  extension  of 
the  separation  shock  from  the  inviscid  flow  region  to  the  flat  test 
plate. 

This  distance,  XP,  was  first  compared  with  the  sum  (AXS  +  XS) , 
for  it  was  believed  that  the  difference  between  these  two- 
dimensional  quantities,  defined  as 

D  =  XP  -  (AXS  +  XS)  (4) 

would  be  such  that 

D-D  (Re/L,  h,  6,  n)  (5) 

where  n  is  the  exponent  of  the  turbulent  boundary- layer  profile. 

The  plotted  results,  shown  in  Figure  15  imply  that  this  difference, 
D,  is  in  fact  independent  of  Re/L  and  h  within  the  range  of  values 
covered  in  this  study.  Further,  since  6  changed  only  fractionally 
throughout  this  experimental  study  and  n  was  nearly  constant  (see 
Appendix  B) ,  it  is  logical  to  assume  that 

D  =  D  (8,  n)  =  constant  =1.0"  (6) 


is  realistic. 

Using  XP  as  a  normalixing  distance,  the  ratios  AXS/XP  and 
XS/XP  may  be  broadly  interpreted  as  the  fraction  of  the  shock 
projection  distance  assigned  to  the  initial  steep  pressure  rise 
process  and  the  fraction  of  the  shock  projection  distance  assigned 
to  the  separated  region,  respectively.  For  the  case  of  the  step, 
Figure  16  shows  the  functional  dependence  of  these  fractional 
distances  on  n/6  for  different  values  of  Re/L.  It  is  not  unreasonable 
to  e;.pect  when  one  also  considers  the  results  of  Werle,  et  al  (Ref. 

(5))  that  as  h/6  increases  without  limit,  that  is,  either  h  becoming 
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large  or  6  becoming  small,  or  both,  AXS/XP  approaches  zero.  For 
the  same  conditions,  XS/XF  approaches  unity.  On  the  other  hand, 
as  h/6  decreases  XS/XP  approaches  zero.  An  example  of  this 
behavior  is  seen  in  the  pressure  distribution  of  Figure  11a  for 
h  =  0.247  inch.  The  limiting  behavior  of  AXS/XP  as  h/<5  approaches 
zero  is  not  known.  It  is  seen  experimentally  that  AXS/XP  increases 
as  h/6  decreases. 

•The  variation  in  the  bluntness  of  the  pressure  distributions 
in  the  free  interaction  region  was  observed  previously  in  Figure  11. 
Since  the  magnitude  of  AXS  is  related  to  this  bluntness  and  the 
bluntness  appears  to  vary  with  Re/L,  one  would  expect,  based  upon 
this  variation,  that  AXS/XP  decreases  with  increasing  Re/L 
(for*  constant  h/6) .  In  the  same  fashion  XS/XP  increases  with 
increasing  Re/L  (for  constant  h/5) .  These  two  facts  are  readily 
verified  in  Figure  15. 

Figure  17  presents  the  geometric  parameters  AXS/h,  XS/'h 
and  XP/h  as  functions  of  Reynolds  number  per  foot.  A  simple 
interpretation  of  the  last  parameter,  XP/h  may  be  made.  This 
parameter  establishes  the  shock  standoff  distance  relative  to  the 
step  and  may  be  interpreted  as  the  cotangent  of  the  turning  angle 
necessary  to  produce  the  oblique  shock  ^ound  in  the  inviscid 
flow  field.  The  proof  of  this  last  statement  is  presented  in 
Figure  18  where  a  comparison  is  made  between  the  experimentally 
measured  plateau  pressures  and  those  computed  using  the  turning 
angle  consistent  with  the  experimental  cotangent,  XP/h.  The  data 
comparison  seems  to  support  the  interpretation  given  XP/h. 

Returning  to  Figure  17,  the  opposite  slopes  exhibited  by 
AXS/h  and  XS/h  are  such  that  XP/h  is  observed  to  be  nearly  constant 
for  a  given  value  of  h — a  result  consistent  with  Equations  (4) 
and  (6) . 

Zukoski  (Ref.  (1))  reports  that  the  wedge  appears  to  be  a 
reasonable  model  for  the  separation  region  when  h  >  1.56.  For 
this  case,  the  cotangent  of  the  wedge  angle  is  4.2  (and  constant) 
for  fully  developed  turbulent  boundary-layer  flow.  Note  in 
Figure  17,  XP/h  is  approximately  6.8  for  a  step  of  height  h  = 

0.603  inch  (n/6  ~  0.2+)  and  6.3  for  a  step  of  height  b  =  1.502 
inches  (h/6  ~  Q.5+).  Hence,  the  experimental  trend  of  XP/h  with 
increasing  h/6  does  not  appear  inconsistent  with  Zukoski’ s  model 
for  h  >  1.56. 


SIDE  FORCE  CALCULATION 


The  induced  control  force  due  to  the  flow  separation  may  be 
evaluated  from  the  upstream  pressure  distributions  as  follows: 
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After  normalizing,  Equation  (7)  may  be  expanded  into  the  following 
form: 
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which  includes  the  two  pressure  integrals,  a  and  8,  of  Figure  14d. 
Introducing  the  shock  standoff  distance,  XP,  and  substituting 
for  the  pressure  integrals,  Equation  (8)  becomes: 


F _  _  P1~P°>  ( XP  \  /  AXS 

T? J~  \h  )  \XP  ' 


(9) 
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where  the  two  integral  expressions,  a  and  8/  have  constant  values 
of  0.400  and  0.969,  respectively. 


i. 


By  a  similar  analysis  based  upon  an  inviscid,  turning 
flow  wedge  model  (the  separation  region  is  replaced  with  a  wedge) , 
the  control  force  may  be  calculated  as: 


Using  the  following  approximations  of  Zukoski  (Ref.  (1) ) 
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Equation  (10)  reduces  to: 


2.1  M 

00 


(13) 


and  is  shown  plotted  in  Figure  19  for  h/S  >  1.5. 

The  control  force  results,  obtained  by  integrating  the 
experimental  pressure  distributions,  are  also  presented  in  Figure 
19  as  a  function  of  h/S.  Note  that  these  experimental  results, 
F/P^h,  increase  with  increasing  Re/L  and  h/6  and  indicate  that 
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for  h/6  >  1.5,  the  dimensionless  force  approaches  the  results 
of  (13)  above.  For  this  approach  to  be  true,  however,  Equation  (9) 
must  reduce  to  Equation  (10) .  Quite  obviously,  two  conditions 
must  be  satisfied:  First,  XP  must  approach  XS  as  h/6  exceeds 
1.5,  and  secondly,  the  ratio  AXS/XP  must  approach  zero  as  h/6 
exceeds  1.5.  These  two  trends  have  already  been  demonstrated  by 
the  data  in  Figure  16.  As  a  result,  the  pressure  coefficient 
term,  (AXS/X?  a  -r  XS  /XP  6),  in  Equation  (9)  also  approaches  1.0 
as  h/6  >  1.5.  Figure  20  simply  verifies  this  fact.  Hence,  Equation 
(9)  reduces  to  Equation  (10)  as  h/6  >  1.5. 

CONCLUSION 

The  pressure  distributions  measured  in  the  separated  region 
ahead  of  steps  were  found  to  be  functions  of  both  Re^  and  h/6 

for  the  turbulent  boundary- layer  separation  case  where  h  <  6 . 

Since  the  induced  side  forces  are  determined  from  these  same 
pressure  distributions,  these  forces  are  also  functions  of  Re^ 

and  h/6.  The  major  result  of  this  study  is  the  definition  of  a 
universal  pressure  distribution  valid  for  two-dimensional  steps. 

Three  geometric  lengths,  AXS,  XS,  and  XP  were  defined.  The 
ratio,  XP/h,  was  interpreted  as  the  cotangent  of  the  wedge  angle 
necessary  to  produce  the  plateau  pressure,  PI.  By  definition, 

XP  locates  the  inviscid  shock  relative  to  the  step  face.  Since 
the  wedge  angle  is  known,  the  shock  angle,  itself,  is  satisfactorily 
computed  using  the  two-dimensional  oblique  shock  relations. 

The  ratios  AXS/XP  and  XS/XP  are  indicative  of  the  fractional 
distances  necessary  for  the  occurrence  of  the  initial  steep  pressure 
rise  and  for  the  extent  of  the  separated  region,  respectively. 

This  last  point  should  perhaps  be  further  qualified.  Zukoski's 
separation  criterion  (Ref.  (1)  as  used  for  this  study  locates  a 
point  on  the  pressure  signature  which  must  be  interpreted  as  a 
scale  change  point,  not  necessarily  a  separation  point  as  the  title 
of  this  criterion  would  suggest.  The  separation  point  per  se 
is  noc  known. 

The  difference  between  XP  and  the  sum,  (AXS  +  XS) ,  while  a 
constant  for  this  study,  is  believed  to  be  functionally  dependent 
upon  the  actual  boundary- layer  thickness,  6,  and  the  turbulent 
v  boundary- layer  profile  exponent,  n. 
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NOTE:  THE  TWO-DIMENSIONAL  WIND  TUNNEL 
NOZZLE  WIDTH  IS  14  INCHES. 


FIG.  I  THE  TWO-DIMENSIONAL  BOUNDARY  LAYER  CHANNEL'S  ADJUSTABLE  NOZZLE 
SHOWING  THE  GLASS-PORTED  SIDE  PLATES  INSTALLED 


INCHES 


FIG.  2  MACH  NUMbER  AND  BOUNDARY  LAYER  THICKNESS  VARIATIONS  V/ITI 
WIND  TUNNEL  SUPPLY  PRESSURE 
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FIG.  5  FLAT  PLATE  SURFACE  PRESSURE  TAP  IDENTIFICATION  NUMBERS  AND  HOLE 
LOCATIONS  RELATIVE  TO  THE  END  PLATES  AND  STEP  MOUNTING  PLANE 


FLAT  PLATE  MOUNTING  INSERT 

A  SECTIONAL  SIDE  VIEW  SHOWING  THE  DETAILS  OF  THE 
STEP  ATTACHED  TO  THE  FLAT  TEST  PLATE 

TEST  STEP  .HEIGHTS  h  =  1.502,  1.400,  1.205,  0.903,  0.603, 

0.247,  AND  0.130  INCHES 

FIG.  6  FORWARD  FACING  STEP  MOUNTING  DETAIL 
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FIG.  8  BOUNDARY  LAYER  CHANNEL  AND  STEP  MODEL  CONFIGURATION 
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STEP  FACE  RELATIVE  TO  Tl 
TWO  SHIFTED  PRESSURE 
DISTRIBUTIONS 


X  DISTANCE 

TYPICAL  REPEATED  PRESSURE  DISTRIBUTIONS  WITH  STEP  FACE  SHIFTED  USING  ALTERNATE  MOUNTING 
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h  =  0.903  INCH  Tq  =  592°R 

P  =  5  ATM  XnM  =  5.50  INCHES 

o  OIL 


h  »  0.603  INCH  Tq  =  592°R 

Pq  =  5  ATM  XQ|L  =  3.50  INCHES 


FIG.  12  OIL  FLOW  SEPARATION  STUDY  (CONTINUED) 
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h  =  1.205  INCHES  Tq  =  592°R 

PQ  =  5  ATM  XQ1L  =  7.70  INCHES 


FIG.  12  OIL  FLOW  SEPARATION  STUDY 
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FIG.  12  OIL  FLO V/  SEPARATION  STUDY  (CONCLUDED) 
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APPENDIX  A- 

FORWARD  FACING  STEP,  TEST  DATA 

This  appendix  lists  the  experimental  test  data  for  the  forward 
facing  step  test  program.  Photographs,  taken  during  the  tests, 
may  be  identified  with  the  corresponding  data  by  means  of  the 
run  number. 

The  reader  should  be  cautioned  that  visible  oil  streaks  (in 
addition  to  the  obvious  crack  in  the  glass)  appear  in  some  of  the 
pictures  cf..  Runs  135  and  193.  The  source  of  this  oil  was  the 
glass  retaining  ring  located  on  the  inside  surface  of  the  .ported 
end  plates.  No  effect  was  observed  on  the  flow  field  as  a  result 
of  this  oil  flow. 

Starting  with  Run  193  and  continuing  through  196,  a  black 
elliptical  shaped  spelt  can  be  observed  on  the  photographs.  This 
spot  was  caused  by  the  condensation  of  high  humidity  wind-tunnel 
room  air  on  the  outside  of  the  wind-tunnel  window.  The  picture 
of  Run  197,  without  the  black  spot,  was  taken  after  the 
condensation  was  evaporated  by  surface  blowing  on  the  outside 
of  the  wind-tunnel  door  window . 
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TABLE  A-l  EXPERIMENTAL  DATA 


RUN  NUMBER 


PARAMETER 

106 

107 

134 

135 

h  (inch) 

1.502 

1.502 

1.205 

1.205 

M 

4.902 

4.675 

4.903 

4.732 

PQ  (psia) 

73.97 

14.04 

73.49 

14.71 

Tq  ( °R) 

592. 

597. 

594. 

597. 

Re/L  x  10“6 
(1/ft) 

4.00 

0.826 

3.96 

0.846 

RUN  NUMBER 


PARAMETER 

150 

151 

166 

167 

h  (inch) 

0.903 

0.903 

0.603 

0.603 

M 

4.911 

4.792 

4.888 

4.764 

PQ  (psia) 

74.04 

15.19 

72.66 

14.47 

Tq  ( °R) 

"  596. 

589. 

604. 

598. 

Re/L  x  10-6 
(1/ft) 

3.96 

0.871 

3.83 

0.819 

RUN  NUMBER 


PARAMETER 

180 

193 

194 

195 

1! 

h  (inch) 

0.247 

1.400 

1.400 

1.400 

M 

4.886 

4.893 

4.881 

4.866 

PQ  (psia) 

72.57 

73.48 

59.26 

43.85 

Tq  ( °R) 

591. 

596. 

594. 

597. 

Re/L  x  10“6 
(1/ft) 

3.97 

3.96 

3.22 

2.38 

A-2 
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RUN  NUMBER 


PARAMETER 

196 

197 

h  (inch) 

1.400 

1.400 

M 

4.834 

4.735 

PQ  (psia) 

28.40 

13.85 

Tq  (°R) 

596. 

602. 

Re/L  x  10“6  (1/ft) 

1.57 

0.784 
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FIG.  A-2,  H  1.205  INCH':0. 
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RUN  180 


FIG.  A-5,  H  =  0.247  INCHES 
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APPENDIX  B 

THE  BOUNDARY-LAYER  DESCRIPTION 

The  description  of  the  boundary  layer  ahead  of  the  step  is 
given  in  this  Appendix.  The  character  of  the  boundary  layer  is 
adequately  described  by  graphs  and  tables  for  the  two  test 
Reynolds  numbers  of  4.0  x  10®  and  0.8  x  10^  per  foot.  .  The 
boundary- layer  data  is  the  result  of  probe  surveys  taken  at 
two  probing  ports,  located  4.8  and  60  inches  downstream,  of  the 
sonic  throat  of  the  NOL  Boundary  Layer  Channel  .{Ref.  (B-l) )  . 

The  probe  tips  were  located  three  inches  ahead  of  the  port 
location.  Therefore,  three  inches  must  be  subtracted  from  the 
48-  and  60-inch  station  values  to  determine  the  true  physical 
location  where  the  probe  data  apply.  The  leading  static  pressure 
tap  (Number  i  in  Fig.  5  of  the  text)  in  the  test  section  was 
located  at  the  76.875-inch  station.  Measurements  made  "at  the 
two  survey  locations  therefore  represent  the  boundary-layer 
conditions  located  19.875,  i.e.,  76.875  less  60.0  plus  3.0,  and 
31.875  inches  upstream  of  this  leading  static  pressure  top 
position. 

The  reasons  for  presenting  these  data  are  two-fold.  First, 
the  character  of  the  boundary  layer  is  identified,  and,  second, 
numerical  data  are  provided  for  those  who  might  wish  to  study 
this  step  separation  problem  by  numerical  methods. 

The  variables  recorded  at  each  probing  port  included  Pitot 
pressure,  equilibrium  temperature  (Ref.  (B-2) )  and  distance. 

The  Boundary  Layer  Channel's  stagnation  pressure  and  temperature 
were  also  recorded.  Data  reduction  followed  procedures  similar 
to  those  described  in  References  (B-3) ,  (B-4)  and  (B-5) .  Since 
the  step  separation  test  was  conducted  in  a  zero-pressure 
gradient  boundary  layer,  Reference  (B-3)  also  serves  as  a 
valuable  source- of  additional  background  material. 
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DELTA 

DELSTAR 

H 

M 

P 

RE/L 

RETHETA 

RHO 

T 

THETA 

THETAE 

THETAH 

U 

y 

Suffix 


FS 

L 

0 

S 

W 


SYMBOLS 

boundary- layer  thickness/  inch 
boundary-layer  displacement  thickness,  inch 
shape  factor,  DELSTAR/THETA 
Mach  number 
pressure,  psia 

unit  Reynolds  number  per  foot 
Reynolds  number  based  oh  momentum 
density,  LBM/FT3 
temperature,  °R 

boundary-lays ^  momentum ‘thickness,  inch 
boundary-layer  energy  thickness,  inch 
boundary- layer  total  enthalpy  thickness ,  inch 
velocity,  ft/sec 

distance  normal  to  plate  surface,  inch 

free  stream 
local 

stagnation 

static 

wall 
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Y 

RHOL 

TSL 

POINT 

INCH 

ML 

PSl/PW 

TOl/TOLFS 

RHOLFS 

TSLFS 

Ul/ULFS 

49 

1.7681 
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1,0060 
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1,0000 
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FIG.  8-2  (U)  TEMPERATURE  DISTRIBUTION:  PO  =  73.  2  PSIA;  TO  =  592  R;  STATION  =  48  INCHES  (U) 


0.1000  0.2000  0.3000  0.4000  0.5000  0.6000  0.7000  0.8000  0.9000 

U/U1NF 

FIG.  B-3  (U)  DIMENSIONLESS  TEMPERATURE-DIMENSIONLESS  VELOCITY  DISTRIBUTION: 
PO  =  73,  2  PSIA;  TO  =  592  °R;  STATION  =  48  INCHES  (U) 
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FIG.  B-5  (U)  LOGARITHMIC  VELOCITY  D 
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FIG.  B-8  (U)  DIMENSIONLESS  TEMPERATURE-VELOCITY  DISTRIBUTION: 

PC  =  15.0  PSIA;  TO  =  588  °R;  STATION  =  48  INCHES  (U) 
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FIG.  B—  1 7  (U)  TEMPERATURE  DISTRIBUTION:  PO  =  14.9  PSIA;  TO  =  586  R;  STATlOh 1  -  60, INCHES  (U> 
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FIG.  B— 18  (U)  DIMENSIONLESS  TEMPERATURE-VELOCITY  DISTRIBUTION: 

PO  =  14. 9  PSIA;  TO  =  586  °R;  STATION  =  60  INCHES  (U) 
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FIG.  B-19  (U)  DIMENSIONLESS  VELOCITY-DISTANCE  DISTRIBUTION: 

PO  =  14.9  PSiA;  70  =  586  °R;  STATION  =  60  INCHES  (U) 


